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"^ ' Abstract 

' ' A complete set of A + 1 mutually unbiased bases (MUBs) exists in Hilbert 

^ I spaces of dimension N = p^, where p is a prime number. They mesh naturally 

^ ' with finite afiine planes of order A, that exist when N = p^. The existence 

O^. of MUBs for other values of A is an open question, and the same is true for 
finite affine planes. I explore the question whether the existence of complete 

'k>( I sets of MUBs is directly related to the existence of finite affine planes. Both 

Vh I questions can be shown to be geometrical questions about a convex polytope, 
but not in any obvious way the same question. 
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1. Introduction 

The acronym MUBs stands for "mutually unbiased bases". By definition, 
two orthonormal bases |ej) and \fi) are mutually unbiased if it is true that 

I (ei I /i) P = constant = — (1) 

for every choice of i,j. Here N is the dimension of the Hilbert space, and 
the second equality gives the value that the constant has to take. Equa- 
tion (1) is reminiscent of the relation between the eigenbases of position 
and momentum — if we have complete knowledge of one observable, we know 
nothing about the other. Bases that are mutually unbiased underlie notions 
such as complementarity. Beginning with Weyl [1], they have been studied 
for a long time [2, 3]. I hope they are of interest also in a conference on 
the foundations of probability and physics. After all, complementarity is one 
of the most striking features of the best probability theory that physics has 
to offer, and MUBs are part of the mathematics of complementarity. The 
importance of MUBs for the foundations of quantum mechanics has been 
stressed, especially, by Brukner and Zeilinger [4]. 

The name MUB is due to Wootters, who pointed out — quite some time 
ago [5] — that if a Hilbert space of dimension A^ contains A^ + 1 bases such 
that every pair of bases is mutually unbiased, then there exists a finite di- 
mensional analogue of the Wigner function, where the role of the classical 
phase space is played by a finite affine plane of order A^. The use of fi- 
nite affine planes in quantum mechanics can be traced back — in the qubit 
case — at least to Feynman [6]; for the most complete treatment so far, see 
Gibbons et al. [7] . Wootters' original motivation had to with quantum state 
tomography. There are other motivations too: Both MUBs and finite affine 
planes play a role in various quantum communication protocols [8] and quan- 
tum cryptography schemes [9, 10] — most spectacularly in a story [11] with a 
happy end [12], known as the Mean King's Problem [13] — but my suspicion 
is that the best is yet to come. Finite geometries are used by people who are 
handling, encrypting and deciphering information classically. With quantum 
information theory in ascendancy it is easy to imagine that, in the long run, 
the appearance of finite affine planes in quantum mechanics will enable us 
to perform a wide range of interesting tricks. A good omen is that some of 
Wootters' key results were rediscovered, and placed in the context of discrete 



mathematics, by Zauner in his Ph D thesis [14]. 

It is known that a Hilbert space of dimension A^ can contain at most 
A^ + 1 MUBs, so we refer to this many MUBs as a complete set. For N = p, 
where p is a prime number, a complete set of MUBs was given by Ivanovic 
[15], and for N = p'' hy Wootters and Fields [16]. A quarter of a century 
later, the question whether a complete set of MUBs exists for other values of 
A^ remains open. Since I have not yet explained what finite affine planes are, 
or what they have to do with MUBs, it may look like an odd fact to mention, 
but it happens that finite affine planes of order A^ exist if A^ = p'^, they do 
not exist when A^ = 6, 10, 14, . . . (an infinite set of values are excluded), and 
it is an open question whether they exist when A^ = 12, 15, 20, . . . (an infinite 
set of values are open). This question goes back to Euler, so after all these 
years it appears justified to say that it is a hard one. Anyway, for discussion 
we have two questions to choose from: 

I: How many MUBs exist in a given Hilbert space? 

II: Given a complete set of MUBs, what can we do with them? 

I choose the first question, since I have nothing to say about the second. I 
do not have an answer to the first either, and in fact I am not sure what kind 
of question it is. We can discuss it as a close packing problem in complex 
projective space [17], or we can try to phrase it as a question of finding 
A^ complex Hadamard matrices related in a particular way — which would 
tie in with what I talked about in last year's Vaxjo meeting [18]. In fact 
the set oi N X N complex Hadamard matrices has quite special properties 
when A^ is prime, or a power of a prime [19]. However, as it stands this 
observation explains nothing — except why I give this talk this year — so I 
choose to explore the obvious question whether the existence of a complete 
set of MUBs is related to the existence of affine planes. The dream is to use 
one of these problems to solve the other! 

I should perhaps say, at the start, that the only conclusion I will arrive 
at is that an affine plane of order A^ exists if and only if it is possible to 
inscribe a regular simplex with A^^ corners in the polytope defined by the 
MUBs (assuming that the latter exist). I would like to think that this sheds 
some light on the question, but I am not so sure. 



2. The Complementarity Polytope 

I will work with density matrices rather than state vectors. An ON basis 
in Hilbert space then becomes a set of A^ projectors Pi = |ej)(ej|. These 
operators sit in the A^^ — 1 dimensional Euclidean space of Hermitian matrices 
of unit trace, with the distance squared between two matrices given by 

D\A, B) = -Tt{A - B)\ (2) 

The factor 1/2 ensures that the projectors in the basis form the A^ corners of 
a regular A^ — 1 dimensional simplex with edges of unit length. Since we will 
encounter several different kinds of simplices, of different dimensionalities, I 
will refer to these simplices as B-simplices. 

We would like to regard our space as a vector space. We can do this by 
placing the origin at the matrix 

The scalar product of two unit trace matrices A and B is then given by 



{A, B) = ] \d\A + B, p,) -D^A- B, p,] 



TyAB - — 

N 



(4) 



The condition for orthogonality of two matrices becomes Ti AB = 1/N, so in 
particular a pair of projectors belonging to two different MUBs are orthogo- 
nal. 

We observe that a complete set of MUBs defines a convex polytope in 
our A^^ — 1 = (A^ — 1)(A^ + 1) dimensional space, which can be described 
very simply: We split our space into A^ + 1 perpendicular A^ — 1 dimensional 
subspaces, and in each subspace we place a regular simplex with A^ corners, 
centered at the origin. Taking the convex cover of all the A^^ + A^ corners we 
obtain a convex polytope that we will call the Complementarity Polytope. 
In terms of matrices, letting the indices i,j run from 1 to A^ and the indices 
/, J run from to A^, the corners define A^(A^ + 1) matrices Pa such that 

TrP,, = TiPl = 1 (5) 



TtPnPij = 0, i^j (6) 



TrPnPjj = 1^ , I^J . (7) 



Equation (5) ensures that all corners sit on a sphere of radius J{N — 1)/2N 
and centered at p^,, equation (6) ensures that they form one regular B-simplex 
for each value of I, and equation (7) ensures that the B-simplices sit in 
N + 1 perpendicular subspaces. Sometimes we will find it convenient to use 
a collective index for the corners, i.e. we will refer to P^j, where u runs 
between 1 and A^^ + A^. 

Reflecting on the definition, we see that it is easy to form a Complemen- 
tarity Polytope, regardless of the value of A^. What is not easy is to ensure 
that its corners are matrices with non-negative spectrum. If they are not, 
they are not density matrices. The question of the existence of a complete 
set of MUBs can now be reformulated: Given a Complementarity Polytope, 
and given the convex body of density matrices, is it possible to arrange them 
in such a way that the former is a subset of the latter? Evidently this is not 
going to be a walkover — the corners of the Complementarity Polytope lie on 
the outsphere of the body of density matrices, and the set of pure quantum 
states there form a 2{N — 1) dimensional subset of this A^^ — 2 dimensional 
sphere. When N > 2 this means that all the corners of the Complementarity 
Polytope must be fitted into a small subset of its outsphere. But we can 
study the Complementarity Polytope regardless of the existence of a com- 
plete set of MUBs, since — as a polytope — it certainly exists in all A^^ — 1 
dimensional vector spaces. This is exactly what we will do, until close to the 
end of the talk. 

The Complementarity Polytope can be visualized easily, at least as long 
as we are not too demanding about how much we want to see. When A^ = 2 it 
is an octahedron. More schematically, we can think of the octahedron as con- 
sisting of three orthogonal one dimensional B-simplices, i.e. as J J J . When 
A^ = 3 we are dealing with four totally orthogonal triangles, A A A A , 
when A^ = 4 we have five totally orthogonal tetrahedra, ^ ^ ^ ^ ^ , 
and so on for any A^. We will be interested in the face structure of our poly- 
tope, and more especially in the following features: When A^ > 2 all edges 
are extremal. There are A^^+^ faces formed by taking the convex cover of 



one corner from each simplex (in the A^ = 3 case this is A A A A , and its 
relatives). These faces are themselves A^ dimensional regular simplices, and 
I call them point faces for a reason that will become clear later on, at least 
to those of you who are familiar with Wootters' phase points [5, 7]. To each 
point face we will associate a point face operator 

point face 

where the index a labels the point faces and the sum runs over all the corners 
in the point face simplex. A more understandable form of the same equation 
is 

Aa-P*= E (^- - P*) ■ (9) 

point face 

It is immediate that 

TiAa = 1 TtAI = N . (10) 

In our vector space the point face operator corresponds to a point on the ray 
from the origin through the center of the point face. If continued in the other 
direction, this ray hits the center of a facet, that is the convex cover of A^ — 1 
corners from each of the B-simplices (in the A^ = 3 case A A A A , and 
its relatives). Hence point faces and facets are in one to one correspondence, 
and they are placed opposite to each other on the surface of the polytope. 
The case A^ = 2 is special since its point faces are also facets, and conversely. 
Note that every point face operator defines a set of parallel hyperplanes 
through the equation 

Tr Ap = constant , (11) 

where p is an Hermitian matrix of unit trace. When p lies in a facet we find 
that Tr Ap = 0, and when p lies in a point face we find that Tr Ap = 1. This 
is a strict proof that the faces that we are looking at are extremal. It also 
proves an interesting conjecture by Galvao [20]. 



3. Finite afRne planes 

Now we come to a seemingly odd observation, namely that the number of 
corners of our polytope — A^^ + N — is equal to the number of lines in a finite 
affine plane of order A^. In the best spirit of Euclid, we give the axioms that 
define affine planes. An affine plane is an ordered pair of two sets, the first of 
which consists of elements a^, called points, and the second of which consists 
consists of subsets l^j of the first, called lines. Two lines whose intersection 
is empty are called parallel. The following axioms hold: 

Al: If Qa and a^ are distinct points, there is a unique line l^^ such that a^ G l^^ 

and a/3 E lu}- 

A2: If tta is a point not contained in the line Z^^, there is a unique line /g- such 

that fla G la and la (^ leu = 0- 

A3: There are at least two points on each line, and there are at least two 

lines. 

To see how this works, think of an ordinary affine plane, and try to think of 
it as two sets, the set of points and the set of lines. Two points determine a 
unique line, while two lines either intersect in a unique point, or else they are 
parallel and do not intersect at all. This is precisely what the axioms say, 
and it is enough to begin to understand how the finite affine planes work. 

We will need only some simple properties. First we define a pencil of 
parallel lines as a maximal set of mutually non-intersecting lines. Then we 
learn, perhaps by looking into some accessible reference [21], that if the 
number of points is finite, then every line contains the same number of points. 
An affine plane of order A^ has A^ points on each line, A^^ points, and N'^ + N 
lines. Each point is on A^ + 1 lines and each pencil of parallel lines contains 
A^ lines. There are altogether A^ + 1 such pencils. 

To construct an affine plane of order A^ we begin by forming a quadratic 
array of A^^ points, representing two pencils of parallel lines. We must then 
somehow find A^ — 1 additional ways of dividing these points into pencils of 
parallel lines, in such a way that a line from one pencil always shares exactly 
one point with a line from any other pencil. For A^ = 2 it is easy: 



1 U W UU WW,. 

1st pencil: 2nd pencil: (12) 



3d pencil: (13) 
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For A^ = 3 we get, after slightly more thinking, 
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(14) 
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2nd pencil: 
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3d pencil: 
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4th pencil: 
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(16) 



(17) 
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We see that an affine plane is a solution to a combinatorial problem. 

As A^ grows it becomes more difficult to ensure that lines from different 
pencils intersect once and once only, so it comes as no surprise to learn that 
it is actually impossible for some values of A^ (such as A^ = 6, 10, 14, . . .). 
Perhaps it is surprising that mathematicians are at a loss to say — even though 
they have had ample time to study the question — whether affine planes exist 
for certain other values of A^ (such as A^ = 12, 15, 20, . . .). The only values 
of A^ for which solutions are known are N = p'^, where p is an arbitrary 
prime number. In fact for prime powers the existence of finite affine planes is 
guaranteed, because for these values of A^ affine planes can be coordinatized 
using finite fields, in just the same way that the ordinary affine plane can be 
coordinatized by real numbers. A field is a set of "numbers" that we can add 
and multiply, in exactly the way we are used to from our experience with 
rational, real, and complex numbers. For any A^ we can do addition and 



subtraction modulo A^, but the existence of an inverse to multiplication is 
non-trivial. It works for N = p, and — after a step resembling that from real 
to complex numbers — also for N = p^, but there are no further possibilities. 
This does not in itself settle the issue of existence of affine planes of order 
N, because there do exist affine planes that are unrelated to fields. This 
happens for instance for A^ = 9, and indeed all known examples occur for A^ 
being a power of a prime [21]. But these examples at least serve to make it 
conceivable that the combinatorics could work out also for, say, A^ = 12. 



4. Painting facets black 

Our question is: Can we associate an affine plane of order A^ with the Comple- 
mentarity Polytope? We can identify the lines with the corners, but exactly 
what collection of A^^ parts of the polytope can serve as points? It turns out 
that what we can do is to select a subset of point faces and declare that they 
are to serve as points. Equivalently — because of the one to one correspon- 
dence between point faces and facets — we can select a subset of A^^ facets, 
paint them black, and declare that the affine plane consists of the black facets 
and all the corners. This is perhaps the most picturesque way to proceed, 
and might be the preferred way of an A^^ — 1 dimensional being — should he 
stoop to consider a finite plane — but for us it is better to work with the point 
faces. The combinatorics will be the same. 

We tackle it by considering the purely geometrical problem of inscribing 
a regular simplex with A^^ corners in the Complementarity Polytope, in such 
a way that each corner of the simplex sits at the center of a point face. We 
call it a D-simplex since it has full dimension. Equivalently, let the corners 
of the D-simplex sit outside the polytope, on points represented by the point 
face operators Aa- We get a regular D-simplex if, restricting the index a to 
range from 1 to A^^, we can choose a set of point face operators so that 

TiA^Afs = N6^fs . (18) 

If so, our selected set of point face operators actually forms an orthonormal 
basis in the A^^ dimensional set of Hermitian matrices, or in other words 
in the Lie algebra of U{N). But can we do it? Recalling the definition 
(8), which tells us that each point face operator represents a choice of one 



corner from each B-simplex, it is not hard to see what the requirement is. 
The iV^ point face operators must be chosen so that for any pair of them, 
the choice of corner must be the same for one B-simplex, and different for 
all the others. This is a combinatorial problem — and the solution of this 
combinatorial problem will turn out to be an afiine plane. 

To begin, let us make choices for the first two B-simplices. Since there are 
only iV^ choices we must use all of them, and we use these choices to organize 
a square array of A^^ points. Using A^ = 3 as our example, the result of this 
preliminary step is 

AAAA AAAA AAAA 

AAAA AAAA AAAA (19) 

AAAA AAAA AAAA 

The understanding is that we must make some definite choice also for the 
A^ — 1 B-simplices that have been left empty here. Let us consider the first 
of these. Interestingly, the condition that we extracted from eq. (18) means 
that we are looking for a Latin square: The problem is to place letters from 
an alphabet of N letters on a square array, in such a way that a given letter 
occurs exactly once in each row and exactly once in each column. There are 
many Latin squares. The difficulties begin when we come to the next step, 
since there are now two Latin squares that force consistency conditions on 
each other. Mathematically, we are facing the problem of finding altogether 
A^ — 1 Latin squares that are mutually orthogonal: By definition a pair of 
Latin squares is said to be orthogonal if every pair of letters, one from each 
Latin square, defines a unique position in the array [21]. 

For the A^ = 3 case, it is enough to find one pair of orthogonal Latin 
squares, and by referring back to the picture of the pencils of parallel lines in 
the affine plane of order A^, we realize that the third and fourth pencil define 
such a pair: 
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AAA AAA 

3d pencil: AAA 4th pencil: AAA (20) 

AAA AAA 

This is not an accident. The A^ + 1 pencils of parallel lines in an affine plane 
of order A^ always define A^ — 1 mutually orthogonal Latin squares — the first 
two pencils are used to organize the array. The converse is also true: An 
affine plane of order A^ exists if and only if one can find A^ — 1 mutually 
orthogonal Latin squares [21]. 

Interestingly, for A^ = 6 one cannot find even one pair of orthogonal Latin 
squares. This was conjectured by Euler and proved 120 years later by Tarry 
[22]. 

To continue our A^ = 3 example, it is clear that our two Latin squares 
can be used to complete the array that we began above, thus: 

AAAA AAAA AAAA 
AAAA AAAA AAAA (21) 

AAAA AAAA AAAA 

This selection of point face operators tells us how to inscribe the D-simplex 
in the Complementarity Polytope. 

The procedure works if and only if an afhne plane of order A^ exists. Any 
affine plane of order A^ will do — we are not restricted to those that can be 
coordinatized with finite fields. If A^ is such that no affine plane of order A^ 
exists, the D-simplex cannot be inscribed in the Complementarity Polytope. 
For values of A^ such that an affine plane exists, we can proceed to reproduce 
all the key formulas in the paper by Gibbons, Hoffman and Wootters [7] . First 
we reindex our selection of point face operators as Aij, where 1 < i,j < N. 
The matrix that represents a corner of the Complementarity Polytope can 
be recovered by summing all the Aij that sit on the line representing that 
corner: 
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^- = ^ E ^^i • (22) 

line LU 

We can also use the fact that our Aij form a basis in the space of matrices to 
define Wootters' Wigner function Wij, which is indeed a function of position 
in the afiine plane: 

W,, = ^TtA.p . (23) 

This equation can easily be solved for p in terms of Wij. Given a Wigner 
function, we can assign a set of numbers p^^ to the lines in such a way that 
Y^Plj = 1, where the sum goes over all lines belonging to a given pencil of 
parallel lines. We do this by summing the Wigner function along the line: 

p^= Y. ^^j = T^^P^P ■ (24) 

line LU 

And SO on [7]. 

Make sure to notice that all of this is true, regardless of the existence 
of the complete set of MUBs. All that is needed is the Complementarity 
Polytope as such, together with an inscribed A^^ — 1 dimensional D-simplex. 



5. MUBs and conclusions 

Seemingly, we have arrived at all the conclusions of Gibbons et al. [7], with- 
out serious work. In fact we have not, because we have tried to settle the 
account without consulting the host. Although the Wigner formulation seems 
to work in happy obliviousness of the question whether the corners are quan- 
tum states or not, there is a catch: For the probability interpretation to 
work, we need Puj > in equation (24). But the only way to ensure this, for 
all positive operators p, is to insist that the P^^ have non-negative spectrum. 
Therefore we must ask whether it is possible to rotate the Complementarity 
Polytope in such a way that it defines a subset of the convex body of density 
matrices. This is our reformulation of the standard problem of defining a 
complete set of MUBs, but it leads to no new insights into it. All I can 
do is to refer you to the literature, where various mathematical approaches 
are employed in clever ways [23, 24, 25, 26, 27] that so far have borne fruit 
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only for A^ = p^. For iV = 6, Zauner's suspicion is that only 3 MUBs exist 
[14], but nobody knows for sure. Wootters has recently reviewed the whole 
question [28]. 

Incidentally, for any polytope in iV^ — 1 dimensions, we can ask if it 
can be oriented so that its corners are pure quantum states. If this can be 
done for a regular D-simplex the result is called a Symmetric Informationally 
Complete POVM [29]. It would be nice to report that the D-simplex that we 
inscribed in the Complementarity Polytope defines a MUB-SIC-POVM — to 
pick an acronym — when the corners of the original polytope become MUBs. 
Unfortunately, at least if we make use of the arguably most elegant choice of 
face point operators [5], this is so only for N = 2 and A^ = 3. 

In conclusion, in the language of the Complementarity Polytope, an affine 
plane of order A^ exists if and only if it is possible to inscribe a regular simplex 
in the Complementarity Polytope, with the iV^ corners of the former situated 
in the middle of the facets of the latter. The problem of finding iV + 1 MUBs 
is equivalent to the problem of rotating the Complementarity Polytope so 
that its corners sit inside the convex body of density matrices. On the face 
of it, these problems are very different. It appears that the finite geometries 
are related to MUBs much like the icing on a cake — it makes the cake more 
delicious, but it has nothing to do with the existence of the cake. This may 
be a superficial view, however. Only the Baker knows. 
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